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Abstract 

There is currently a great deal of interest in the theoretical and practical possibility 
of cloaking objects from the observation by electromagnetic waves. The basic idea 
of these invisibility devices [HI El H3] , [18] is to use anisotropic transformation media 
whose permittivity and permeability e Xu , fi Xu , are obtained from the ones, e Xu , [i Xv , 
of isotropic media, by singular transformations of coordinates. 
In this paper we study electromagnetic cloaking in the time-domain using the for- 
malism of time-dependent scattering theory [23] . This formalism provides us with 
a rigorous method to analyze the propagation of electromagnetic wave packets 
with finite energy in transformation media. In particular, it allows us to settle 
in an unambiguous way the mathematical problems posed by the singularities of 
the inverse of the permittivity and the permeability of the transformation media 
on the boundary of the cloaked objects. Von Neumann's theory of self-adjoint 
extensions of symmetric operators plays an important role on this issue. We 
write Maxwell's equations in Schrodinger form with the electromagnetic propa- 
gator playing the role of the Hamiltonian. We prove that the electromagnetic 
propagator outside of the cloaked objects is essentially self-adjoint. This means 
that it has only one self-adjoint extension, A^, and that this self-adjoint extension 
generates the only possible unitary time evolution, with constant energy, for finite 
energy electromagnetic waves, propagating outside of the cloaked objects. 

*PACS classification scheme 2006: 41.20.Jb, 02.30.Tb,02.30.Zz, 02.60.Lj. 
* Research partially supported by CONACYT under Project P42553F. 
■^Fellow Sistema Nacional de Investigadores. 
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Moreover, Aq is unitarily equivalent to the electromagnetic propagator in the 
medium e Xp ,fi Xu . Using this fact, and since the coordinate transformation is the 
identity outside of a ball, we prove that the scattering operator is the identity 
This implies that for any incoming finite-energy electromagnetic wave packet the 
outgoing wave packet is precisely the same. In other words, it is not possible 
to detect the cloaked objects in any scattering experiment where a finite-energy 
wave packet is sent towards the cloaked objects, since the outgoing wave packet 
that is measured after interaction is the same as the incoming one. Our results 
give a rigorous proof that the construction of [HI EJ [13] , [18] cloaks passive and 
active devices from observation by electromagnetic waves. Actually, the cloaking 
outside is independent of what is inside the cloaked objects. 
As is well known, self-adjoint extensions can be understood in terms of bound- 
ary conditions. Actually, for the electromagnetic fields in the domain of Aq the 
component tangential to the exterior of the boundary of the cloaked objects of 
both, the electric and the magnetic field have to be zero. This boundary condi- 
tion is self-adjoint in our case because the permittivity and the permeability are 
degenerate on the boundary of the cloaked objects. 

Furthermore, we prove cloaking for general anisotropic materials. In particular, 
our results prove that it is possible to cloak objects inside general crystals. 

1 Introduction 

There is currently a great deal of interest in the theoretical and practical possibility of 
cloaking objects from the observation by electromagnetic fields. The basic idea of these in- 
visibility devices [HI M, E32] , [IB] is to use anisotropic transformation media whose permittivity 
and permeability, e Xu ,fi x ", are obtained from the ones, e Xv ,i1q u , of isotropic media, by sin- 
gular transformations of coordinates. The singularities lie on the boundary of the objects to 
be cloaked. Here the material interpretation is taken. Namely, the e Xu ,fi X " and the e Xu ,fi Xu , 
represent the components in flat Cartesian space of the permittivity and the permeability of 
physical media with different material properties. It appears that with existing technology it 
is possible to construct media as described above using artificially structured metamaterials. 
In [HI [9] a proof of cloaking was given for the conductivity equation -i.e., in the case of zero 
frequency- from detection by measurement of the Dirichlet to Neumann map that relates 
the value of the electric potential on the boundary to its normal derivative. The papers [13] 
and [18J consider electromagnetic waves in the geometrical optics approximation, i.e. for 
large frequencies. In [23] a experimental verification of cloaking is presented and [3] and 
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[5] give a numerical simulation. A rigorous prof of cloaking has already been given by [7] 
where fixed frequency waves were studied, i.e., in the frequency domain. They consider a 
class of finite energy solutions to Maxwell's equations in a bounded set, O, that contains the 
cloaked object on its interior, and they prove cloaking, at any frequency, with respect to the 
measurement of the Cauchy data of these solutions on the boundary of O. We give further 
comments on this paper below. For other results on this problem see [25] and [T5]. In [TB] 
cloaking of elastic waves is considered, and the history of invisibility is discussed. 

In this paper we study electromagnetic cloaking in the time-domain using the formalism of 
time-dependent scattering theory [23]. This formalism provides us with a rigorous method to 
analyze the propagation of electromagnetic wave packets with finite energy in transformation 
media. In particular, it allows us to settle in an unambiguous way the mathematical problems 
posed by the singularities of the inverse of the permittivity and the permeability of the 
transformation media on the boundary of the cloaked objects. Von Neumann's theory of self- 
adjoint extensions of symmetric operators plays an important role on this issue. We write 
Maxwell's equations in Schrodinger form with the electromagnetic propagator playing the 
role of the Hamiltonian. We prove that the electromagnetic propagator outside of the cloaked 
objects is essentially self-adjoint. This means that it has only one self-adjoint extension, 
An, and that this self-adjoint extension generates the only possible unitary time evolution, 
with constant energy, for finite energy electromagnetic waves propagating outside of the 
cloaked objects. Moreover, Aq is unitarily equivalent to the electromagnetic propagator in 
the medium Eq u , [Iq v . Using this fact, and since the coordinate transformation is the identity 
outside of a ball, we prove that the scattering operator is the identity. This implies that for 
any incoming finite-energy electromagnetic wave packet the outgoing wave packet is precisely 
the same. In other words, it is not possible to detect the cloaked objects in any scattering 
experiment where a finite-energy wave packet is sent towards the cloaked objects, since the 
outgoing wave packet that is measured after interaction is the same as the incoming one. 
Our results give a rigorous proof that the construction of [HI [9], [13] , [18] cloaks passive and 
active devices from observation by electromagnetic waves. Actually, the cloaking outside is 
independent of what is inside the cloaked objects. 

As is well known, self-adjoint extensions can be understood in terms of boundary condi- 
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tions. Actually, for the electromagnetic fields in the domain of Aq the component tangential 
to the exterior of the boundary of the cloaked objects of both, the electric and the mag- 
netic field have to be zero. This boundary condition is self-adjoint in our case because the 
permittivity and the permeability are degenerate on the boundary of the cloaked objects. 

Furthermore, we prove cloaking for general anisotropic materials. In particular, our 
results prove that it is possible to cloak objects inside general crystals. 

Even though, as mentioned above, the cloaking is independent of the cloaked objects, 
and in particular, the cloaking outside is not affected by the presence of passive and/or 
active devices inside the cloaked objects, we discuss the dynamics of electromagnetic waves 
inside the cloaked objects for completeness, since it helps to understand the above mentioned 
independence of cloaking from the properties of the cloaked objects. 

We prove that every self-adjoint extension of the electromagnetic propagator in a trans- 
formation medium is the direct sum of the unique self-adjoint extension in the exterior of 
the cloaked objects, Aq, with some self-adjoint extension of the electromagnetic propagator 
in the interior of the cloaked objects. Each of these self-adjoint extensions corresponds to a 
possible unitary time evolution for finite energy electromagnetic waves. As is well known, 
the fact that time evolution is unitary assures us that energy is conserved. This results 
implies that the electromagnetic waves inside and outside of the cloaked objects completely 
decouple from each other. Actually, the electromagnetic waves inside the cloaked objects 
are not allowed to leave them, and viceversa, the electromagnetic waves outside can not go 
inside. 

In terms of boundary conditions, this means that transmission conditions that link the 
electromagnetic fields inside and outside the cloaked objects are not allowed, since they do 
not correspond to self-adjoint extensions of the electromagnetic propagator, and then, they 
do not lead to a unitary dynamics that conserves energy. Furthermore, choosing a particular 
self-adjoint extension of the electromagnetic propagator of the cloaked objects amounts to 
choosing some boundary condition on the inside of the boundary of the cloaked objects. In 
other words, any possible unitary dynamics implies the existence of some boundary condition 
on the inside of the boundary of the cloaked objects. 
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The fact that there is a large class of self-adjoint extensions -or boundary conditions- 
that can be taken inside the cloaked objects could be useful in order to enhance cloaking 
in practice, where one has to consider approximate transformation media as well as in the 
analysis of the stability of cloaking. 

Actually, we consider a slightly more general construction than the one of [HI El [13] , [18] 
since we allow for a finite number of cloaked objects. 

In [7j a very general construction for cloaking is introduced. In the case of Maxwell's 
equations all their constructions are made within the context of the permittivity and the 
permeability tensor densities being conformal to each other, i.e., multiples of each other by a 
positive scalar function. In particular, all isotropic media are included in this category. They 
mention that both for mathematical and practical reasons it would be very interesting to 
understand cloaking for general anisotropic materials in the absence of this assumption. In 
this paper we actually solve this problem, since we prove cloaking for all general anisotropic 
materials. In particular, our results prove that it is possible to cloak objects inside general 
crystals. 

Note, moreover, that [7J also considers the cases of the Helmholtz equation. We do not 
discuss this problems here. 

Furthermore, remark that the existing theorems in the uniqueness of inverse scattering 
do not apply under the present conditions. 

In [7J cloaking is proven with respect to the Cauchy data at any fixed frequency given 
on a surface that encloses the cloaked object. In the case where the permittivity and the 
permeability are bounded above and below it is well known that the Cauchy data at a fixed 
frequency is equivalent to the scattering matrix at the same frequency. See for example 
[T7] and [2B]- This equivalence is, however, not proven in the case where the permittivity 
and the permeability are degenerate on the boundary of the objects. In fact, it is perhaps 
even not true for general degenerate media that are not transformation media since in this 
case it is possible that there are finite energy electromagnetic waves that are absorbed by 
the boundary of the objects as t — > ±oo. If this is true, the equivalence will not hold 
since the Cauchy data in a surface that encloses the objects will not contain information 
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on the waves that are asymptotically absorbed by the boundary of the objects. It is a 
problem of independent interest to see if this actually happens or not for general degenerate 
permittivities and permeabilities. For an example of scattering by a bounded obstacle with 
a singular boundary and Neumann boundary condition, where this happens see [TO]. For a 
similar situation in the scattering of electromagnetic waves by a Schwarzschild black-hole see 
[2] • Note that in our approach we directly consider the scattering operator that is measured 
in scattering experiments. 

In the analysis of Maxwell's equations with permittivity and permeability that are inde- 
pendent of frequency the dispersion of the medium is not taken into account. This means 
that cloaking will hold for electromagnetic wave packets with a narrow enough range of 
frequencies, such that this assumption is valid. 

The paper is organized as follows. In Section 2 we prove our results in electromagnetic 
cloaking. In Section 3 we consider the propagation of electromagnetic waves in the interior 
of the cloaked objects. In Section 4 we formulate cloaking as a boundary value problem 
outside of the cloaked objects for the Maxwell equations at a fixed frequency, following 
our analysis of the self-adjoint extensions of the electromagnetic propagator. In particular, 
we give the appropriate boundary condition on the outside of the boundary of the cloaked 
objects. Finally, in Section 5 we prove cloaking of infinite cylinders. This is of interest 
since this is the case considered in the experimental verification in [21] and in the numerical 
simulations of [1] and [5]. Of course, [21] only consider a slice of the cylinder. In Sections 3 
and 4 we give further comments on the results of [7]. 

Addendum 

After the previous version of this paper was posted in the arXiv we published the paper 
[2"5] where we generalized the results of this paper on spherical cloaks to the case of high- 
order cloaks, and where we also discussed cloaking in the frequency domain. Moreover, in 
our paper [30] we identified the cloaking boundary condition that has to be satisfied in the 
inside of the boundary of the cloaked objects, in the case where the permittivity and the 
permeability are bounded above and below inside the cloaked objects. 



6 



2 Electromagnetic Cloaking 



Let us consider Maxwell's equations, 



VxE = -| B , vxh=|d, (2.1) 



VB = 0,V-D = 0, (2.2) 
in a domain, O C I 3 , as follows, 

n := M 3 \ uf =1 ^, Kj n = Q,j^l (2.3) 

where = 1, 2, • • • , N, are the objects to be cloaked. We assume that each Kj is a ball 
with center Cj and radius Oj, i.e., 

A} = {x G M 3 : |x - Cj\ < dj} , j = 1, 2, • • • , N. (2.4) 
The cloaked objects are denoted by 

K := uf =l K y 

We designate the Cartesian coordinates of x by x x , A = 1, 2, 3 and by E\, H\, B x , D x , X = 
1, 2, 3, respectively, the components of E, H, B, and D. As usual, we denote by e Xu and fi Xu , 
respectively, the permittivity and the permeability. We have that, 

D x = e Xu E u , B x = yL Xu H v , (2.5) 

where we use the standard convention of summing over repeated lower and upper indices. 

We consider now a transformation from Q := R 3 \ {c 1; c 2 , ■ ■ • , Cjy} onto Q that was first 
used to obtain cloaking for the conductivity equation, i.e. at zero frequency, by [HI E] and 
then by [18J for cloaking electromagnetic waves (for a related result in two dimensions using 
conformal mappings see 
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For any y 6 K 3 we denote, y := y/|y|. Let y x , X = 1,2,3, designate the cartesian 
coordinates of y G fi . Take bj > dj,j = 1, 2, • • • , N. Then, for < |y — Cj\ < bj, we define, 

x = x(y) = /(y) := Cj + -c,\ + a,- J y^. (2.6) 

Note that this transformation blows up the point Cj onto and that it sends the punctu- 
ated ball Be jib j) := {y G 1R 3 : < |y — Cj| < bj} onto the spherical shell, a,j < |x — Cj\ < bj. 
We assume that, 



B Cj (bj) n B Cl (bi) = Q,j?l,l<j,l< N. (2.7) 

For y G M 3 \U^ =1 5c we define the transformation to be the identity, x = x(y) = /(y) := 
y. Our transformation is a bijection from Q onto Q. By y = y(x) := / _1 (x) we designate 
the inverse transformation. We denote the elements of the Jacobian matrix by Ay, 



A X> ■= 7^37- (2-8) 



dx 



dy x ' ' 

Note that the A\, G C 1 U} \ \jf =1 dB Cj (bj)\ We designate by A% the elements of the 
Jacobian of the inverse bijection, y = y(x) := / _1 (x), 

Af^^e&fa^dBcthj). (2.9) 

The papers [SJ E] and [TS] considered the case where N — 1, Ci = 0. 

We take here the so called material interpretation and we consider our transformation 
as a bijection between two different spaces, Qq and Q. However, our transformation can be 
considered, as well, as a change of coordinates in fl . Of course, these two point of view 
are mathematically equivalent. This means, in particular, that under our transformation the 
Maxwell equations in Q and in Q will have the same invariance that they have under change 
of coordinates in three-space. See, for example, [2T]. Let us denote by A the determinant of 
the Jacobian matrix (12.81) . Then, 

2 

, forO < \y-Cj\ < bj. (2.10) 



bj-aj ( -VHy-Cjl + aj 



b j V |y - 9 
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This result is easily obtained rotating into a coordinate system such that, y — Cj = (|y — 
Cj |,0,0) |25J. Fbryefio\UjL 1 Sc i (6 J -),A = l. 

Let us denote by E , H , B , D , e Xu , {j, Xv , respectively, the electric and magnetic fields, 
the magnetic induction, the electric displacement, and the permittivity and permeability of 
Qq. Eq 1/ ,[Iq, are positive, Hermitian matrices that are constant in Q . 

The electric field is a covariant vector that transforms as, 

E x (x) = AX{y)E 0tX ,(y). (2.11) 

The magnetic field H is a covariant pseudo- vector, but as we only consider space transfor- 
mations with positive determinant, it also transforms as in (12. lip . The magnetic induction 
B and the electric displacement D are contravariant vector densities of weight one that 
transform as 

5 A (x) = (A(y))- 1 ^(y)5 A '(y), (2.12) 

with the same transformation for D. The permittivity and permeability are contravariant 
tensor densities of weight one that transform as, 

e *"(x) = (A( y ))-^,(v) ^,(v)4 v (y), (2-13) 

with the same transformation for yr v . The Maxwell equations (\2.1\ 12.21) are the same in 
both spaces Q and fl . Let us denote by e Xv , n\ u , e 0Xu , /j,q Xu , respectively, the inverses of the 
corresponding permittivity and permeability. They are covariant tensor densities of weight 
minus one that transform as, 



eju,(x) = A(y)^'(yX(y)5 0AV (y), ^(x) = A(y)ylf(y)<(y) tev ,(y). (2.14) 
Note that 

det e Xu = A' 1 det e Xu , det fi Xu = A' 1 det (2.15) 

det e Xv = A det e 0Xu , det fj, Xu = A det fi 0Xu . (2.16) 
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We now introduce the Hilbert spaces of electric and magnetic fields with finite energy. 
The E , H , B , D , were defined in fl , but since M. 3 \ Q = {cj}^ =1 is of measure zero, we 
can consider them as defined in M 3 , what we do below. 

We denote by Hoe the Hilbert space of all measurable, square integrable, C 3 — valued 
functions defined on M 3 with the scalar product, 

(e«",E« 2 >) := / (2-17) 

V /OE J R3 

We similarly define the Hilbert space,Ho# , of all measurable, square integrable, C 3 — valued 
functions defined on M. 3 with the scalar product, 

(H«Hf) := / HS^dy'. (2.18) 

V / OH J R 3 

The Hilbert space of finite energy fields in M 3 is the direct sum 

H :=Hoe®Hoh- (2.19) 

Moreover, we designate by Hne the Hilbert space of all measurable, C 3 — valued functions 
defined on Q that are square integrable with the weight e Xu , with the scalar product, 

( E (D )E (2))^ ;= f E ^e Xv Wd^. (2.20) 
Jn 

Finally, we denote by H-cih the Hilbert space of all measurable, C 3 — valued functions 
defined on Q that are square integrable with the weight n Xu , with the scalar product, 

(H«,H( 2 ))^:= / H^^Wd* 3 . (2.21) 
Jn 

The Hilbert space of finite energy fields in Q is the direct sum 

H a :=H QE ®Hn H - (2.22) 

We now write the Maxwell's equations (12.1 1) in Schrodinger form. We first consider the 
case of M. 3 . We denote by Eq and u, , respectively, the matrices with entries eo\ u and fJ,o\ v . 
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Recall that (V x E) = s Xl/p -£-E p where s Xvp is the permutation contravariant pseudo- 
density of weight —1 (see section 6 of chapter II of [2T], where a different notation is used). 
By ao we denote the following formal differential operator, 



/E \ ./ £„VxH„ \ 
(Hoj =l ( -/i V x E J • ^ 

Here, as usual, we denote, EqV x H := £q\ u (V x Hq) u , and /xoV x E = /ioA^(V x Eq) i/ . 



Then, equations (12.11) are equivalent to, 
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Let us denote by Cj(IR 3 ) the set of all C 6 — valued continuously differentiable functions 
on M. 3 that have compact support. Then, ao with domain Cq(IR 3 ) is a symmetric operator in 
Ho, i.e., a C ag. Moreover, it is essentially self-adjoint in TC , i.e., it has only one self-adjoint 
extension, that we denote by A . Its domain is given by, 



Eo \ ( E 



D(A ) = {( H Q ) : G M H ' ^ f ' ( 2 - 25 ) 



Eq \ / Eq \ / Eq 



and, 

^(5:j=°°(5:j'(Hj efl (^ < 2 - 26 » 

where the derivatives are taken in distribution sense. These results follow easily from the 
fact that -via the Fourier transform- ao is unitarily equivalent to multiplication by a matrix 
valued function that is symmetric with respect to the scalar product of H.q. Moreover, it 
follows from explicit computation that the only eigenvalue of Aq is zero, that it has infinite 
multiplicity, and that, 



ft ± := (kernel Ao) 1 = j ( ^° J G H : ^o^ou = 0, ^o^ou = J . (2.27) 

Furthermore, A has no singular-continuous spectrum and its absolutely-continuous spec- 
trum is R. See, for example, [271 128] . 

Taking any 

' GK oi n £>(4,) (2.28) 



E 
H 
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we obtain a finite energy solution to the Maxwell equations (\2.1\ 12. 2p as follows 



(2.29) 



This is the unique finite energy solution with initial value at t — given by (12.281) . Note 
that as e~ ltA °Ho± C Ho± equations (12. 2p are satisfied for all times if they are satisfied at 



Let us now consider the case of fl We denote by e and /i, respectively, the matrices with 
entries E\ v and fi\ u . 

We now define the following formal differential operator, 



Let us denote by Cq(^2) the set of all C 6 — valued continuously differentiable functions on 
f2 that have compact support. Then, with domain Cl(fl) is a symmetric operator in Hq. 
To construct a unitary dynamics that preserves energy we have to analyse the self-adjoint 
extensions of oq. 

We denote by Ue the following unitary operator from TL QE onto Hqe, 



t = 0. 




(2.30) 



Equations (12. ip are equivalent to, 




(f^Eo) A (x) :=A^E oy (y), 



(2.31) 



and by Uh the unitary operator from Hqh on to TCqh, 



(^H ) A (x) :=A^H 0X/ (y). 



(2.32) 



Then, 



U := U E © Ui 



(2.33) 



is a unitary operator from Hq onto Hq. 
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We denote by a 00 the restriction of a to Cj(f2 ). The operator a 00 is essentially 
self-adjoint and its only self-adjoint extension is Aq. This follows from the essential self- 
adjointness of ao and from the fact that any function in Cq(1R 3 ) can be approximated in the 
graph norm of ao by functions in Cj(fio). To prove this take any continuously differentiable 
real-valued function, 0, defined on R such that, <fi(y) = 0, \y\ < 1 and 4>(y) = 1, \y\ > 2. 
Then, for any 

( Ho ) e c ^ 

we have that, 

n^Ny-c.l) (h ) GC o( fi o) 

and moreover, 




where by s-lim we designate the strong limit in Ti,Q. 

As aoo is essentially self-adjoint, it follows from the invariance of Maxwell equations that 
an is essentially self-adjoint, and that its unique self-adjoint extension, that we denote by 
An, satisfies 

A n = UA U*. (2.34) 
For the proof of these facts see [29]. Hence, we have the following theorem. 

THEOREM 2.1. The operator is essentially self-adjoint, and its unique self-adjoint 
extension, A n , satisfies {2.31$ . 

The unitary equivalence given by (I2.34p implies that Aq has the same spectral prop- 
erties that Aq. Namely, it has no singular-continuous spectrum, the absolutely-continuous 
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spectrum is R and the only eigenvalue is zero and it has infinite multiplicity. Moreover, 




(2.35) 



(2.36) 




(2.37) 



This is the unique finite energy solution with initial value at t — given by (12.361) . Note 
that as e~ ltAn Ti.Q± C H.n± equations (12.21) are satisfied for all times if they are satisfied at 
t — 0. We can consider more general solutions by considering the scale of spaces associated 
with An, but we do not go into this direction here. 

The facts that is essentially self-adjoint and that its unique self-adjoint extension Aq is 
unitarily equivalent to the propagator A Q of the homogeneous medium are strong statements. 
They mean that the only possible unitary dynamics in Q that preserves energy is given by 
(12.371) and that this dynamics is unitarily equivalent to the free dynamics in R 3 given by 
(12.291) . In fact, dQ acts like a horizon for electromagnetic waves propagating in Q in the 
sense that the dynamics is uniquely defined without any need to consider the cloaked objects 
K = L)j =1 Kj. As we will prove below this implies electromagnetic cloaking for all frequencies 
in the strong sense that the scattering operator is the identity. 



Since D{A n ) = UD{A Q ), for any (E, H) T e D{A n ) there is a (E , H ) T G D{A ) such 



that 




(2.38) 



Then, it follows from (I23T1 [27321 [2T33D that 



E x n = 0, H x n = 0, in dK. 



(2.39) 
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where dK + denotes the outside of the boundary of the cloaked objects, K, and n is the 
normal vector to dK + , if (E ,H ) are, for example, bounded near dK + . That is to say, for 
electromagnetic fields in the domain of the tangential components of both, the electric 
and the magnetic field vanish in the exterior of the boundary of the cloaked objects. This is a 
self-adjoint boundary condition because the permittivity and the permeability are degenerate 
on dK + . 

Let xn be the characteristic function of fl, i.e., x^( x ) = 1, x G O, Xfi( x ) — 0, x G IR 3 \ Q. 
We define, 



J( £ j j (x) := xn(x) I h q )(*)• (2-40) 



By (EHEnnira), 



|e Al, (x)|<C, |/i A!y (x)| < C, X6S1. 

Then, J is a bounded operator from Hq into 7Yq. 
The wave operators are defined as follows, 

W± = s- lim e itAn Je- itAo P 0±} (2.41) 

where Po_i_ denotes the projector onto Ho±- 

Let us designate by W 1,2 (IR 3 ) the Sobolev space of C 6 valued functions. We denote by / 
the identity operator on TLq. Then, 

LEMMA 2.2. 

W± = UP 0± . (2.42) 

Proof: Denote, 

W(t):= e itAn Je~ UAo P ±. 

By (E31, for any <p G Tio 

W{t)<p = iP{t) + UP 0± <p, (2.43) 
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with 

iP(t) := Ue itAo (U*J - I) e- itAo P 0± <p. (2.44) 

Let Br denote the ball of center zero and radius R in M. 3 . Since for \y\ > R, with R 
large enough, our transformation, x = /(y), is the identity, x = y, and in consequence, 
^v(y) = ^ or |y| — P j we nave that, 



(U*J-I) = {U*J-I) XB . (2.45) 



It follows that, 



with, 



s- lim ip(t) = Us- lim e itA °${t) (2.46) 

t— >±oo t^ioo 



;= (fT J - /) XBR e- itA °P 0±¥ >. (2.47) 



We have that, 



||tf(*)|| Wo < JxB R e ° p o±<P H + XB R e ° p o^ 



< C 



XB R e- ttA °Po±V ■ (2.48) 

rto 



Then, as (A + i)~ 1 P ± is bounded from H into W 1,2 (R 3 ) [27] [28], it follows from the 
Rellich local compactness theorem that 



01 



is a compact operator in Hq. Suppose that if G D(Aq) fl 7^o±- Then, 



s- lim x Bn e ~ UAop o±V = s- I™ X^JA) + P^e~ UA ° {A, + i)^ 

t— >±oo ti f^+oo -^K 



and whence, by (12.481) . 



s- lim ■Qit) = 0, 

t— >±oo 



(2.49) 



(2.50) 
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and it follows that in this case, 



s- lim ifj(t) = 0. (2.51) 

t— >±oo 

By continuity this is also true for ip G 'Ho±. 
Then, (12^421) follows from (g33D and (12311 . 

□ 

The scattering operator is defined as 

S:=W*iy_. (2.52) 

COROLLARY 2.3. 

5 = P 0± . (2.53) 
Proof: This is immediate from (12.421) because [/* U — I. 

□ 

Let us denote by S± the restriction of S to 7io±- <Sj_ is the physically relevant scattering 
operator that acts in the Hilbert space 7i _i_ of finite energy fields that satisfy equations (12. 2p . 
We designate by I± the identity operator on Hq±. We have that, 

COROLLARY 2.4. 

S ± = I ± . (2.54) 

Proof: This follows from Corollary 2.3. 

□ 

The fact that S± is the identity operator on TCo± means that there is perfect cloaking for 
all frequencies. Suppose that for very negative times we are given an incoming wave packet 
e~ ltA °(p-, with y?„ G TCo±- Then, for large positive times the outgoing wave packet is given 
by e~ ltA °ip + with ip + = S_iLp_. But, as S± = I, we have that <p + = ip_ and then, 

e - itAo ip_ = e- itAo ip + . 
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Since the incoming and the outgoing wave packets are the same there is no way to detect 
the cloaked objects K from scattering experiments performed in Q. 

In this paper we considered transformation media obtained from a singular transformation 
that blows up a finite number of points, by simplicity, and since this is the situation in the 
applications. Suppose that we have a transformation that is singular in a set of points 
that we call M and denote now Qq := IR 3 \ M. What we really used in the proofs is that 
W 1,2 (R 8 ) = Wj' 2 (n ) where Wj' 2 (O ) denotes the completion of Cg°(n ) in the norm of 
W 1,2 (IR 3 ). We also assumed that e$ v , iXq v are constant. What was actually needed is that a® 
is essentially self-adjoint. All our results hold under this more general conditions provided 
that in (12.411 |2~4"2"1) and (12.531) we replace Pq± by the projector onto the absolutely-continuous 
subspace of A Q and that we assume that D(A ) n Ho ac C W 1,2 (IR 3 ), where we have denoted 
the absolutely-continuous subspace of A Q by Ho ac - Moreover, S± has to be defined as the 
restriction of S to Ho ac and in (I2.54p I± has to be the identity operator on Ho ac - Note that 
under these general assumptions A could have non-zero eigenvalues and singular-continuous 
spectrum. 

For example, W li2 (R 3 ) = Wj' 2 (n ) if M has zero Sobolev one capacity [Tl [TTl H2|. 
Moreover, assume that the permittivity and the permeability tensor densities Eq , ^ are 
bounded below and above. Under this condition ao is essentially self-adjoint. Furthermore, 
let us denote by Ho the Hilbert space of finite energy solutions defined as in (12.191) but 
with Eq u = /Zq M = d~ XfJl . Let A ,Ho± be, respectively, the electromagnetic propagator in Ho 
and the orthogonal complement of its kernel. We have that Ho and Ho are the same set 
of functions with equivalent norms. Furthermore, D(Aq) = D(Aq), kernel Aq = kernel Ao. 
Moreover, (E , H ) T G Ho± if and only if E = Eo^o, Ho = /ioH for some (E , H ) G Ho±- 
As [23 [28] D(A ) n H ± C W 1 ' 2 (M 3 ) we have that D(A ) n H ± C W 1 - 2 (M 3 ) if e , fx are 
bounded operators on W 1,2 (IR 3 ) and this is true if the derivatives ^-^o, ~^~^o are bounded 
operators on Ho for p = 1, 2, 3. Note, furthermore, that Ho ac C Hq±. 
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3 Electromagnetic Waves Inside the Cloaked Objects 



Let us now consider the propagation of electromagnetic waves in the cloaked objects. For 
this purpose we assume that in each Kj the permittivity and the permeability are given 
by Ej,fj,j U , with inverses ej\ u ,fij\ u and where £j,jij are the matrices with entries Ej\ v ,^j\v 
Furthermore, we assume that < Ej, fij u < C, x G Kj and that for any compact set Q con- 
tained in the interior of Kj there is a positive constant Cq such that det ej u > Cq, det \xf > 
Cq, xgQ. In other words, we only allow for possible singularities of Sj, fij on the boundary 

Of Kj. 

We designate by TljE the Hilbert space of all measurable, C 3 — valued functions defined 
on Kj that are square integrable with the weight Ej V , with the scalar product, 

(E«Ef) j£ :=^£<^^x'. (3.1) 

Similarly, we denote by TCjH the Hilbert space of all measurable, C 3 — valued functions 
defined on Kj that are square integrable with the weight fij U , with the scalar product, 

^J^tflfli*. (3.2) 
The Hilbert space of finite energy fields in Kj is the direct sum 

TCj := 7~LjE © 7~tjH-> (3-3) 
and the Hilbert space in the cloaked objects K is the direct sum, 

Hk '■— ®f =1 Hj. 

The complete Hilbert space of finite energy fields including the cloaked objects is, 

H:=H Q ®H K - (3.4) 
We now write (12.11) as a Schrodinger equation in each Kj as before. We define the 
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following formal differential operator, 



a, { %)=i( I. (3.5) 



H i / V x E i 



Equation (12. lft in f^,- is equivalent to 



Let us denote by Co(i^j) the set of all C 6 — valued continuously differentiate functions on 
Kj that have compact support in the interior of Kj, that we denote by Kj := Kj\ dKj. 
Then, a,j with domain Cq(Kj) is a symmetric operator in TCj. We denote, 



a := a n ® j=1 CLj, (3.7) 

with domain, 

D{a) := { ( H„ ) ®^ ( ^ ) G ®^ =1 C ° ( ^ } } • (3 ' 8) 

The operator a is symmetric in Ti. The possible unitary dynamics that preserve energy for 
the whole system including the cloaked objects K are given by the self-adjoint extensions of 
a. Let us denote a the closure of a, with similar notation for an, aj,j — ,N. Then, 

a = A n ®f =1 a], 

where we used the fact that as an is essentially self-adjoint, an = An- The adjoint of a is 
given by, 



D{a*) 



Ef2 \ A 



'- [ la |eH: 



and 



for 



En 



H, 



En 



Ej 
H, 



He 



Ej 
H, 



(3.9) 

(3.10) 
(3.11) 
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Let us denote by Kq± : = kernel(i T«q), JCj± := kernel(z T a *j) the deficiency subspaces of 
an and aj, j = 1, • • • ,N. Since an is essentially self-adjoint K,q± = {0}. Let K,± := ffi^/Cji 
be the deficiency subspaces of a K := ®f =1 aj. Suppose that JC± have the same dimension. 
Then, it follows from Corollary 1 in page 141 of [22] that there is a one-to-one correspondence 
between self-adjoint extensions of ok and unitary maps from /C + into /C_. If V is such a 
unitary, then the corresponding self-adjoint extension A KV is given by, 



Hence, since /Cq± = {0} and a = © 7Tk there is a one-to-one correspondence between self- 
adjoint extensions of a and unitary maps, V, from K, + into /C_. The self-adjoint extension 
A v corresponding to V is given by, 

A v = A n © A KV . 

Thus, we have proven the following theorem. 

THEOREM 3.1. Every self-adjoint extension, A, of a is the direct sum of Aq and of some 
self-adjoint extension, Ak of ok, i-e., 



This theorem tells us that the cloaked objects K and the exterior Vt are completely 
decoupled and that we are free to choose any boundary condition inside the cloaked objects K 
that makes ax self-adjoint without disturbing the cloaking effect in Q. Boundary conditions 
that make Ak self-adjoint are well known. See for example, [19], [20] , [H] and [6]. 

It follows from explicit computation that zero is an eigenvalue of every Ak with infinite 
multiplicity and that, 



D(A KV ) = + ip + + Vip+ : (f E D(aW), cp + G /C+} 



and 



A = A Q © A K - 



(3.12) 



TCk± '■= (kernel A K ) C 




E 
H 



) 








V 



dx x 







VkH u = 



(3.13) 
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where by e% (x) := e^(x) for x e K,-, and /x^(x) := /i^(x) for x <E i^-, j = 1, 2, • • ■ , N. It 
follows that zero is an eigenvalue of A with infinite multiplicity and that, 

H± : = (kernel A) L = H n ±®H K ±. (3.14) 

For any ip = ip Q © ip K e H± H D(A), 



e - itA <p = e~ itAn cpn © e~ itAK ip K (3.15) 

is the unique solution of Maxwell's equations (12. 11 12.21) with finite energy that is equal to <p 
at t = 0. This shows once again that the dynamics in Q and in K are completely decoupled. 
If at t = the electromagnetic fields are zero in Q, they remain equal to zero for all times, 
and viceversa. Actually, electromagnetic waves inside the cloaked objects are not allowed to 
leave them, and viceversa, electromagnetic waves outside can not go inside. This implies, in 
particular, that the presence of active devices inside the cloaked objects has no effect on the 
cloaking outside. In terms of boundary conditions, this means that transmission conditions 
that link the electromagnetic fields inside and outside the cloaked objects are not allowed. 
Furthermore, choosing a particular self-adjoint extension of the electromagnetic propagator 
of the cloaked objects amounts to choosing some boundary condition on the inside of the 
boundary of the cloaked objects. In other words, any possible unitary dynamics implies 
the existence of some boundary condition on the inside of the boundary of the cloaked 
objects. The particular boundary condition that nature will take depends on the specific 
properties of the metamaterial used to build the transformation media as well us on the 
properties of the media inside the cloaked objects. Note that this does not mean that we 
have to put any physical surface, a lining, on the surface of the cloaked object to enforce any 
particular boundary condition on the inside, since as we already mentioned this plays no role 
in the cloaking outside. It would be, however, of theoretical interest to see what the interior 
boundary condition turns out to be for specific cloaked objects and metamaterials. These 
results apply to the exact transformation media that we consider on this paper. However, 
the fact that there is a large class of self-adjoint extensions -or boundary conditions- that can 
be taken inside the cloaked objects could be useful in order to enhance cloaking in practice, 
where one has to consider approximate transformation media as well as in the analysis of 
the stability of cloaking. 



22 



The fact that for the single coating there has to be boundary conditions on the inside 
of OK has already been observed by [TJ. In Definition 4.1 of [7] a definition of finite energy 
solutions is given. Furthermore, is proven in Theorem 6.1 that in the case of the single 
coating -where the permittivity and the permeability are bounded above and below inside 
the cloaked object- the tangential components of the electric and the magnetic field of these 
solutions have to vanish in the inside of the boundary of the cloaked object. Note that in 
this case in order to have a self-adjoint extension of the electromagnetic propagator inside 
the cloaked object we are only allowed to require that either the tangential component of E 
or the tangential component of H vanishes, but not both. 

These boundary conditions are called hidden boundary conditions in [7] where also the 
case of the Helmholtz equation is considered. In the case of Maxwell's equations they propose 
two solutions to this issue. One of them is a lining, i.e., a physical material on the boundary of 
the cloaked object that enforces a particular boundary condition, for example, they propose 
a lining by a perfect electric conductor. Note that this raises now the question of what is the 
boundary condition between the lining and the cloaking metamaterial. In fact, we face the 
same problem as before, since we can always consider that the lining is part of the cloaked 
objects, and then, the question of what is the appropriate boundary condition remains. The 
second proposal of [7j is a double coating that corresponds to surrounding both the inner 
and the outersurface of the cloaked objects with appropriately matched metamaterials. As 
our permittivities and permeabilities inside K are allowed to vanish as they approach OK 
the double coating fits in our formalism. 

In Theorem 5.1 of [7] cloaking is proven for all frequencies and active devices, with the 
double coating, with respect to the Cauchy data of the finite energy solutions that they define 
in Definition 4.1. 

Remark that there is no real contradiction between our results and the ones of [7J. Our 
results imply that there is always a hidden boundary condition on the inside of the boundary 
of the cloaked objects, that is imposed upon us by the fundamental principle of the con- 
servation of the energy of the electromagnetic waves, that implies that time evolution has 
to be given by a unitary group generated by a self-adjoint extension of the electromagnetic 
propagator, and this amounts to a boundary condition at the inside of the boundary of 
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the cloaked objects. Note that we do not exclude here the possibility that in some cases 
the electromagnetic propagator of the cloaked objects could be essentially self-adjoint, and 
in this situation the dynamics inside the cloaked objects will be uniquely defined. In this 
case the hidden boundary condition will be uniquely determined by the boundary conditions 
satisfied by the functions in the domain of the unique self-adjoint realization of the elec- 
tromagnetic propagator in the cloaked objects. Note, however, that we have proven that 
for exact transformation media the cloaking outside is actually independent of the cloaked 
objects. 

4 Cloaking as a Boundary Value Problem 

It is a question of independent interest to consider cloaking as a boundary value problem for 
the Maxwell' system at a fixed frequency 

V x E = iXB, V x H = -iAD, A ^ 0, (4.1) 
VB = 0,V-D = 0. (4.2) 

As we have already shown, cloaking is independent of the cloaked object, and this means 
that we only have to consider these equation in Q. The main question now is to decide what 
is an appropriate class of solutions with locally finite energy. Our analysis of the self-adjoint 
extensions of the electromagnetic propagator shows that we have to take solutions that are 
locally in the domain of Aq, that is to say that they are given by f!2.38j) 

(SM2) 

with (E , H ) T locally in the domain of Aq, i.e., (E , H ) T are in the domain of A when 
multiplied by any function in C^°(IR 3 ). It follows from (I2.39P that the solutions with locally 
finite energy have to satisfy the boundary condition, 

E x n = 0, H x n = 0, in dK + , 
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where dK + is the outside of the boundary of the cloaked object. This is the only self- 
adjoint boundary condition on dK + . Note that we define in the same way solutions with 
(locally) finite energy in a bounded subset of Q. In [7] a different definition of solutions with 
(locally) finite energy is given in Definition 4.1. 

5 Cloaking an Infinite Cylinder 

We discuss now the case of an infinite cylinder. For simplicity we consider one cylinder 
centered at zero and with its axis the vertical line L := (0, 0, x 3 ), x 3 G R. Then, 

K:={x= (x\x 2 , x 3 ) G M 3 : v 7 ^ 1 ! 2 + \x 2 \ 2 < a, x 3 G m} , fi := M 3 \ K. (5.3) 
The set Q is now given by, 

tt = R 3 \ L. (5.4) 

Let us denote by x := (a;, 1 ,^ 2 ) the vectors in the x 1 — x 2 plane and x := x/|x|. The 
transformation (12.61) is replaced by 

x = x(y) = /( y ):={^=<7fel +a >i' (5.5) 

[ ^3 — 2/3, 

for < |y| < b and with b > a. This transformation blows up the line L onto OK and it 
sends Kb\L onto K^\K where 

K b := {y = (y\y 2 ,y 3 ) e M 3 : vW+W < b,y 3 e m} . 

For y G M. 3 \ K b we define the transformation to be the identity, x = y. 

The Hilbert spaces of finite energy electromagnetic fields, the unitary operator U and a , 
A , an, dR, a , are defined as in Section 2. 

THEOREM 5.1. The operator an is essentially self-adjoint, and its unique self-adjoint 
extension, An, satisfies 

An = UA U*. (5.6) 
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Proof: The theorem is proven as Theorem 2.1 observing that W 1 ' 2 (R 2 ) = W M (R 2 \ 0) since 
{0} has zero Sobolev one capacity in IR 2 [TJ [TTJ [T2] . 

□ 

We now consider the wave and the scattering operators. For simplicity we assume below 
that Eq u = £5 Xu ,fiQ U = fi5 Xu . The wave operators are defined as in ( 12.4=1 jl but now the 
operator J is defined as follows, 



J ( V H JJ (X) : = Xn(x)0(x) j(x), 

where <fi is continuous and it satisfies 0(x) = (|x| — a), a < |x| < a + 5 and 0(x) = 1 for 
|x| > a + 25, for some S > 0. 

LEMMA 5.2. 

W± = E7P ±. (5.7) 

Proof: The lemma is proven as in the proof of Lemma 2.2, but in (12.451 12.461 12.471 12.481) we 
have to replace xr d , by xr n where, Cr := {y G I 3 : lyl < R} for R large enough. Now we 
can not prove (12.491 12.501 12.511) by compactness arguments because K is unbounded. Instead 
we use propagation estimates for A . The following results are well known. See for example 
[21 [271 1211 [31] where the general anisotropic case is considered. For any ip G 7Y ±, 



e - itAo (p = — / e ik - y (e- iuJ+(k)t P + (k)0(k) + e- iuJ - {k)t P.(k)<f(k)) d 3 k (5.8) 

where (p is the Fourier transform of cp, to±(k) = ±\k\c with c := (e/i)~ 1//2 , and P±(k) are 
projectors on R 3 that are infinitely differentiable for k e M 3 \0. Suppose that (p G C^°(1R 3 \L) 
and let O be a bounded open set such that O C M 3 \ L and support (p C O. Denote 

Then by the (non) stationary phase Theorem (see the Corollary to Theorem XI. 14 of [23] ), 
for any n — 1, 2, • • ■ there is a constant C n such that 

|(e"^V)(y)| <C n (i + |y| + |t|r\ ±|^0. (5.9) 

26 



We write, 



with 



and 



Xc fl e-^V = 0i + 02 (5.10) 
>i := X i±y/{cm d)Xc R e- UAo <p (5.11) 



z--=X {±yme 6)Xc R e- itA °ip. (5.12) 



by O 



s- lim 0x = 0. (5.13) 

t— >±oo 

Note, furthermore, that there is an e > such that |fe| > e for any k £ O. Then, for any 
±^ G O, |y| > c|t|e. It follows that there is a T such that 

02 = 0, for|t| > T. (5.14) 

By QSMEMEIM 

s- lim xc R e~ itAo ¥ = 0, (5.15) 

and (12.501 12.511) follow. Note that P _i_ is not needed because <p G 7^ ±- By continuity 
this is true for all if G H ±. Then, (ETJ) follows from (g3Sj I23TI) . 

□ 

The scattering operator is defined as in (12.521) . 
COROLLARY 5.3. 

5 = P 0± . (5.16) 
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Proof: This is immediate from (15 .7p because U* U = I. 

□ 

Let us denote by S± the restriction of S to Hq±. S± is the physically relevant scattering 
operator that acts in the Hilbert space 7^ _i_ of finite energy fields that satisfy equations (12. 2p . 
We designate by I± the identity operator on TCqj_. We have that, 

COROLLARY 5.4. 

S± = I±. (5.17) 

Proof: This follows from Corollary 5.3. 

□ 

Again, the fact that Sj_ is the identity operator on TLqi_ means that there is cloaking for 
all frequencies. 

In Theorem 7.1 of [7J cloaking is proven for all frequencies with respect to the Cauchy 
data of the finite energy solutions that they define in Definition 4.1 and furthermore, in 
Theorem 8.2, they prove cloaking for all frequencies with the SHS boundary condition with 
respect to the Cauchy data of the finite energy solutions that they define in Definition 8.1. 

Theorem 3.1 remains true in the case of the cylinder. The proof is the same. Furthermore, 
all the remarks about finite energy solutions, and cloaking and that we made in Sections 
2, 3, are true in the case of a cylinder. We do not repeat them here. Moreover, equations 
(12.381) hold. However, since now the transformation (15.51) only acts on the plane orthogonal 
to the axis of the cylinder equations (I2.39|) has to be replaced by 



Exx = 0, Hxx = 0, indK + , (5.18) 
where E:= (E 1 ,E 2 ),U-= (#i,# 2 )- 

As in Section 4 we define solutions to (14.11 14.21) with locally finite energy as solutions 
that are locally in the domain of Aq, that is to say that they are given by (I2.38P 



E )=U E ° 
H U 1 H 
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with (E , H ) T locally in the domain of A Q , i.e., (E , H ) T are in the domain of A when 
multiplied by any function in C^°(IR 3 ). It follows from (I5.18P that the solutions with locally 
finite energy have to satisfy the boundary condition, 

Exx = 0,Hxx = 0, in9K + . (5.19) 

Note that (15.191) is the SHS boundary condition considered in [7j. We have proven here 
that (15.191) is the only self-adjoint boundary condition on dK + allowed by energy conserva- 
tion. 
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